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The quantum-mechanical problem of N fermions with 5-function interaction in a one-dimensional 
potential well of finite depth is solved. It is shown that there exists exact wave function of Bethe- 
ansatz form in the case that a single particle tunnels outside of the well. The Bethe-ansatz like 
secular equations for the spectrum are obtained. 
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I. INTRODUCTION 



One-dimensional quantum-mechanical systems of electrons have been a subject of great importance in recent years 
because they serve as tractable limits of quasi-one-dimensional systems as well as models for the conductive 
properties of the one-dimensional quantum wire . A model Hamiltonian describing the system is a A-body problem 
■ with delta-function interaction. For convenience, let fi and the mass of the electron be one. The Hamiltonian of the 
electron gas in one dimension then reads 



> . N q2 N 

<N I H = -J2 jra +c S ( x i~ x j) 



dxj 

= 1 1 3>i=l 



where c stands for the interaction strength. The quantum-mechanical problem of the above Hamiltonian were solved 
in [4-8] with periodic boundary conditions, and in || with the boundary condition for a potential well of infinite depth. 
Preceding a very narrow and very deep attractive potential at one side of the well Ref. ^ considered another kind of 
boundary condition. In the discussion of ||, a phase shift between the incident and reflected waves, which depends 
on a phenomenological parameter, was introduced. The problem of many particles with (S-function interaction in a 
potential well of finite depth is more complicated due to the fact that the wave function outside of the well is nonzero 
and should be taken into account. In lid] one of the author solved the problem of A bosons with delta- function 
interaction in a potential well of finite depth. The more important physical problem should be N fermions in a 
potential well of finite depth. So we will discuss this problem in the present paper. The Hamiltonian of spin-1/2 
A^-particle system interacting by a (5-function in the presence of a square well is 



N o2 N N 



£ : h = -T,q^ + E y (^) + 2c E s(xi- Xj ), (i.i) 

i=l i i=l i>j—l 



where 



V( Xi 



i . L 

M < j 



V a 2 \ Xi \ > 1 
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in which L is the width of the well. When Vq = 0, the Hamiltonian is known to be exactly solvable by the Bethe 
ansatz in the case of periodic boundary condition. When Vb ^ 0, the Hamiltonian is not invariant under translation, 
and the total momentum of the system is not conserved. Howeve r, t he system is still invariant under the action of 
the permutation group Sm- We will diagonalize the Hamiltonian (|Ll|) by a generalized Bethe ansatz which involves 
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all non-diffraction scattering waves. When the system is considered as a canonical ensemble in which almost all the 
particles are confined in the well, except for one particle which is able to tunnel out of the well, we will show that the 
diagonalization can be exactly realized and the consistent Bethe ansatz equations obtained. 

In section || we consider the two-particle system to illustrate the idea of the Bethe ansatz jnj. In section |l| 
we analyze the generic system of N fcrmions and solve the wave function of Bethe ansatz form. In section |j"v| , we 
diagonalize the S'-matrix and obtain secular equations for the spectrum. In section [v], We give some discussions and 
remarks. 



II. TWO-PARTICLE STATES 



Before exploring the generalized Bethe ansatz for a system of N fermions, we discuss the case of two fermions. The 
quantum-mechanical problem of one particle in a potential well of finite depth is well known. Because there is no 
external field, the spin degrees of freedom a = +, — arc degenerate. The wave functions ip a (x) are 

f ^e™, x < -k. 

il> a (x) = \ A^e lkx + A-e- lkx , \x\ < f. (2.1) 
[ A^e~ KX , x > §. 

where k > and k 2 — k 2 + Vq = E is the energy of the system. The amplitudes A's are related via the boundary 
conditions at x = ±-j, which say that the wave function and its derivative are continuous at the ends of the well. 

One can easily find that the spectrum is determined by ( K — ^f ) 2 = e~ 2lkL . 

k — ik 

Because of the <5-function interaction between the two particles, we should divide the {x\, X2) plane into two regions, 
namely, the region C{id) = {(£1,0:2) \x\ < X2} and the region C(a\) = {(x\,X2) \x 2 < x\}. The definition of a\ and 
other terminology will be given in the next section. In the square well x\,x 2 E [—■§■,■§■], the wave function takes the 
form of plane waves of all the non-diffraction scatterings. We suppose that in the region C(id), the wave function 
reads 

r a i a2 (x u x 2 )= [p> felKl+ ^ 2X2) ^a 1 a 2 (a,/3)-e^^ + ' 3fc ^ 1 ' J B a2ai (a,/3)]. (2.2) 

a,0=+,- 

The amplitudes A(a, 0), B(a 1 f3) are so chosen that in the region C(<7i), the wave function reads 

C a >i'*2)= [e i(afclBI+ ^« ) Bai« a («,/^-e i ( afcia!a+ ' 3fe,a!1 ^« !1 a 1 («,/3)]. (2.3) 

a,P=+- 

Obviously, the requirement of antisymmetry for the fermion wave function 4 , ^ ia2 ( Xl ^ X2 ) = ~ V^ai i x 2> x i) ^ s satisfied 
by(U)andQ. 

When the interaction strength is not zero, the wave functions is continuous while its derivative has a finite difference 
at the (5-function singularity 

d d 

fe - g^)«a 2 -<a 2 )( x l, x 2)\ Xl = X2 = 2c^{x X ,X 2 )\ Xl = X2 . (2.4) 
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Substituting (2.2) and (2.3) into (2.4), we obtain a relation between the amplitudes A and B 

B aia2 (a,P) = Sal^(a, f3)A a , ia , 2 (a, 0) (2.5) 

, Q a' x a' 2l a\ (o^h - /3k 2 )S2 1 1 Sa 2 +2icS2 2 1 S2 1 2 . , , , ,, . 

where baia 2 {pt, p) — : — is tne two-body scattering matrix. 

aki — pki + 2ic 

When one of the fcrmions moves out of the potential well, there is an additional restriction that the wave func- 
tion should vanish at infinity. Thus the wave function on the left side of the well takes the form V'aia (^l: x %) — 
E/3=+,- [e^e^ k —A L aia2 (+,/?)- e K "^e i ^"B^ ai {+,l3)] for x x < -f < x 2 < f, and ^ ia2 \x u x 2 ) = 

T l p=+,-[ eK2X2et0klXlB a 1 a 2 (+, P) ~ e KlX2 e if3k2Xl A% 2ai (+,l3)} for x 2 < -f < x x < -§. The former is the antisym- 
metric counterpart of the later and vice versa. Similarly, the wave function on the right side of the well should 
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be < Q2 (zi,z 2 ) = Ef3=+,-[^ K2X2 ^ klXl A^ ia2 (-,p) - e-'^eV^B*^-,®] for -| < Xl < % < x 2 , and 
<J^i,x 2 ) = Ep= + ,-le-^e^ k ^BK a2 (-,{3) - e -^ e ^»A« (1 (- /})] for -f < x 2 < f < a*. 

The amplitudes A fl:i (a, /3), B R ' L (a, (3) are related to the amplitudes /3), B(a.(3) via the boundary conditions 
at x = ±4, namely 

^a l03 (— 2^2) = VVa^"^,^), 
^-^a 1 o 2 ( a; l> a: 2)|xi=-L/2 = ^j-V'aia 2 (a;i,a:2)U 1= _L/2 



and 



Thus we get 



^-V'a ia2 (a:i,a;2)U 2 =L/2 = ^-^a ia2 (xi, x 2 )| X2= i/ 2 . 



iklL Ki-iki _ A aia2 (-,/3) 



( + ,/?)' 



- lk2 L K 2~ik 2 _ B a2ai (/3,-) 



k 2 -H/c 2 S a2ai (/3,+)' 



Kl+iki B a201 (-,/?)' 



and 



c - lk2L K2-ik 2 = A aia2 (f3,+) 
K 2 +ik 2 A aia2 (p, — ) ' 

where /3 = +, — . We observe that the dependence of A' R ' i , i?- R ' i on A, B are the same as what occurred in the one 
particle situation. 



The equations (2.6-2.9) define the secular equations for the spectrum of quasi momentum ( charge ) k\, k 2 . This is 



due to the fact that 

S _1 (fci, k 2 )S{k u -fe) = -S-\-k 2 , kJSfaM (2-10) 
The secular equation for the spectrum is 

S- 1 ^, +)S(/3, -)A(fi, +) = / 2 2 A(/3, +), /? = +,-. (2.11) 
where /j = — — '^-e~ lkiL , j _ ^2. After diagonalizing the scattering operators S and S~ x , one can obtain the 



Bethe ansatz equations for k\,k 2 . This will be discussed in section IV 



III. THE CASE OF N FERMIONS 

We first introduce some terminology which is helpful for avoiding the ambiguities which have appeared in some 
previous literature. In a Euclidean space IR W with Cartesian coordinates x = (xi, x 2 , ■ ■ ■ , xn), the set of hyperplanes 
— Xj — 0}| i,j — 1, 2, • • • N} partition R into finitely many regions. We use the convention that the scalar 

product of two ./V-dimensional vectors is written as (x\y) = ^2^ = i(x)i(y)i. The above-mentioned hyperplanes are 
Weyl reflection hyperplanes (x\atij) = where the roots of the Lie algebra A^-i- We denote a 
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Weyl reflection hyperplane as P a := {x \ (x\a) = 0}, in which a is a root of Am-i- The connected components of 
R N \ {IPq} are called Weyl Chambers (l^] of the Lie algebra Aff—i- We denote the Weyl group of Ajy_i as Wa, the 
basic elements of which are defined by <Ji : (x%, X2, ■ ■ ■ , Xi, Xj+i, • • • , xn ) * (x%, X2, • • • , Xi+\, Xi, • • • , xn )• The Weyl 
group of Bff, denoted by Wb, has one more basic element <jn besides the previously defined = 1, 2, • • • N — 1). 
The definition of cr^v is given by ctjv : {xx, ■ ■ ■ , xjv-i, xjv) | — > (zi, ■ • ■ ,#iv-ij — a?iv) - An element of the Weyl group 
obviously maps one Weyl chamber onto another and all the chambers can be obtained from any given chamber via 
the actions of the whole Weyl group. Thus Weyl chambers can be specified by the elements of Weyl group. 

Now we consider the case of N fermions. Obviously, the Hamiltonian (1.1) is invariant under the action of the 
permutation group Sn, but it is not invariant under tra nsla tion. Thus the total momentum of the system is not 
conserved. The Schrodinger equation of the Hamiltonian (1.1) on the domain R^^Pq} becomes 



N 

E 



d 2 

~^+V(x i ) 



1p a (x) = Elp a (x), 



where a := (ai, 02, • • • , ajy) and eij G {+, — }. Solutions of the equation are plane waves. Considering that the total 
momentum is not conserved, we adopt the following Bethc |11| ansatz form 



4>a{x) 



E 

<7<EW B 



A a (a,T)e l ^ x \xeC(T) 



(3.1) 



where ak stands for the image of a given k := (k\, k%, • • • , fc^r) by a mapping a G Wb and the coefficients A(a, r) 
are functionals on Wb ® Wa- We emphasize that the sum runs over the Weyl group of the Lie algebra Bm but the 
wave function is defined on various Weyl chambers corresponding to the Weyl group of the Lie algebra An-i- This 
is different from the situation of periodic boundary condition. 

For a fermionic system, the wave function is supposed to be anti-symmetric under any permutation of both coor- 
dinates and spin states, i.e. 



(3.2) 



Here (<Jiip) a is well defined by tp ai a(o' i 1 x). Therefore both sides of ( |3.2[ ) can be written out by using (3.1). Furthermore, 
using the evident identity ( ak\a ^ 1 x) = (<7iak\x) and the rearrangement theorem of group theory, we obtain the 
following consequence from (|3.2|) 



A a (<J,<7iT) = -A aia ((Ji<T,T) 



(3.3) 



The (5-function term in the Hamiltonian ( |l . 1[ ) contributes a boundary condition at hyperplane P a (a is a root of 
Lie algebra An-i), namely a discontinuity of the derivative of wave function along the normal of Weyl hyperplane: 



lim [a ■ Vip a (x( a ) + ea) 



lim [ipafaa) + ea ) ~ i>a(%(a) - ea )\ = °> 
c->0+ 



(3.4) 



where xi^ G P a an d V := y^,-_i e^djdxi). 
Substituting (3.1) into (3.4), we find that 

i[(ak)i - (ak) i+1 ][A a (a,aiT) - A a (a i a,a l T) - A a (a,r) + A a (aia, r)] 



2c[A a (<r,r) + A a (ai<r,T]. 



By making use of the relation (^3), we can obtain the following relations from ( p.q ) 

A a {a l a,T) = S l aa ,(ak)A a >(a,T). 



(3.5) 



(3.6) 



St al (cTk) 



c8 a ,a> ~ i[{ak) l - (ak) i+1 }P a 



c — i[(ak)i — (ak) 



i+U 



(3.7) 



where a' = aia and P a . a > stands for the matrix elements of the spinor representation of permutation gr oup . 

The relation ( |3.3| ) provides for the coefficients A a relation between different Weyl chambers. Eq.(3.6) provides 
a connection between those coefficients which are related via any element of Weyl group Wa in the same Weyl 
chamber. Although the 2 N N\ coefficients are determined by (3.6) only up to 2 N arbitrary constants, we will see 
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in next section that the secular equation for the spectrum is determined uniquely. The basic elements of the Weyl 
group Wa obey of = 1 and <Jiai + i<j; L = erj_|_i<Xj<Tj_|_i as identities. These id entit ies must involve some relations 
i.e. A a (afa,r) — A a (a,r) and A a (<7i<7i + i(Ti<T,T) — ^(a-j+icrjC^+icr, r). Using ( p^ ) repeatedly, one can obtain the 
following relations: 

S i (a i ak)S i (ak)=I, 

S' l (a l+1 a l( 7k)S l+1 (<T t ak)S l {ak) = S <+ Vi^m^)SVi+i^)S i+1 MO> (3.8) 

where we have adopted the conventions S = matrix (S a b), S l = S ® I, S 1+1 = I S ( J is a 2 x 2 unit matrix). 
These relations are consistency conditions for the S-matrix. The second relation is called a Yang-Baxter equation. 
The concrete S-matrix in ([3.7]) verifies these relations. 



IV. DIAGONALIZATION OF THE S-MATRIX AND THE SECULAR EQUATION FOR THE 

SPECTRUM 



Because the antisymmetric property provides a relation for wave function between different Weyl chambers, we 
only need to discuss the problem in one chamber. For simplicity we consider C(id) i.e. the region x\ < x 2 < ■ ■ ■ < 
xjv- Owning to the potential well terms in Hamiltonian, we must consider the wave function in the three different 

L L L L , L L 

regions x% < — — < x 2 < ■ ■ ■ < x^ <—,——< %i < x 2 < ■ ■ ■ < xjy < — and — — < x\ < x 2 < ■ ■ ■ < Xff-i < — < xn 

Lj Li Li Li Li Li 

respectively. The cases of more than one particle being outside of the interval [— — , — ] are not important since those 

L L 
regions are not next to the region — — < xi < x 2 < ■ ■ ■ < xn < —. 

We first make some nota- 

tion conventions: W' := {02, 03, • • • , <Jn-i, <tn } and VV" := {01, a 2 , ■ ■ ■ , on-2, cn-iOn^n-i} are two subgroups 
of the Weyl group of Bn; two particular cycles are q'j : (x\, ■ ■ • , Xj, ■ ■ ■ , xn) ► (xj, x\, • • • , Xj-%, Xj+%, • • • , xn) and 
q'j : (x\, • ■ ■ , Xj, ■ ■ ■ , xn) !— > (x%, • • • , Xj-i, Xj+i, • • • , xn, Xj). Then we can write (3^) into a sum of pairs and get the 
following expressions ( because there will be no confusion here, we adopt the notation if> a instead of ipff ): 

f-S T ■ L L 

(I) m the region — — < x\ < x 2 < ■ ■ ■ < xn < — , 

N 

j = l cr'eW 

N 

= E E [M°"lh id y kiXN +A a ^No"q';,id)e- l ^ XN y (a "^ k \ x " ) (4.1) 

j=l cr"eW" 

here and the following x' :— (0, x 2 , • ■ • , x^), x" := (xi,x 2 , ■ ■ • , xjv-i, 0) and <j\ := <J\(Ji ■ ■ ■ om-i&n^n-i • 1 • 02 01; 

r , • L L 

(II) In the region x\ < — — <x 2 <---<xn<— , 

N 

E A L a {o'q' v id)e K ^e^'^ k ^ (4.2) 

j=l o-'eW' 

here and the following Kj > 0, Kj + k^ = U 2 ; 

/. ' /. 

(III) In the region — — < X\ < x 2 < ■ ■ ■ < xn-i < — < xn, 

N 

^( X )=J2 E Af(a"q' J !,id)e- K ^e i ^" q 'i k ^"l (4.3) 

j=l cr"eW" 

Because of the finite depth of potential well, both the wave function and its derivative are continuous at the ends of 
the well. At the left end, we have 
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9a; x 



xi=-L/2 



d 
dxi 



4>a 



x 1 = -L/2 



Thus we get A h a (o' q'^id) = A a (a'q' J ,id)e-^- K ^ L ^ 2 + A a (a 1 a / q' J ,id)eS' lk ^+ K ^ L ^ and 

Kj + ikj A a (a'q' j: id) 
Similarly, the conditions at the right end read 



(4.4) 



(4.5) 



1p a \x a =L/2 = ^a\x 2 =L/2 



They give rise to A^(a"q'J ,id) = A a (a"q'J, id)e^ +K ^ L / 2 + A a (a N a"q'j, id)e^ ik ' +K ^ L l 2 and 

«j + » fc j _ . . MoN<j"g"M) ikjL 



A a (a"q",id) 



(4.6) 



(4.7) 



By making use of ( f4.l|) and (3.6), we obtain the matrix relations 

A(aq'j,id) = Si{a 2 ■ ■ ■ a N -iaq'j'k)S 2 (cr 3 ■ ■ ■ a N -iaq"k) 
■ ■ ■ S N - 2 (crN-icrqjk)S N -i(aq'-k) ■ A(<jq" , id) 

A{diaq' 3 ,id) = Sx(a 2 ■ ■ ■ a N ^id N aq'-k)S 2 {(Jz ■ ■ ■ a N -ia N aq"k) 
■ ■ ■ S N - 2 {<T N -ia N <Tq''k)S N ^i{a N (Tq''k) ■ A(a N aq" ,id). 



(4.8) 



(4.9) 



Here a nd in th e following equ ations th e spin labels are omitted, hence they should be understood as matrix equations. 
From (4.5,4.7) together with (4.8,4.9) we obtain an eigenequation in the spin space: 



where 



and 



Tf%A(a N a^f, id) = f 2 A{a N aq'J, id) 
Tj = iS-" -1 • • • S^S^ N ■ ■ ■ S^^ 1 Tj — S-^^ 1 ■ ■ ■ S^S-* 1 * ■ ■ ■ S^ +1 



S 1 ' = Si(a i+ i ■ ■ ■ <j N -i<rq"k) 



S Jl = Si(a i+1 ■ ■ ■ a N ^x(7 N aq'!k) 



[(ck)j — (ak)i]P — ic 
{ok)j — (ak)i + ic 

[-((Tfc)j - (cTk)i]P-ic 
—(ak)j — {uk)i + ic 



(4.10) 
(4.11) 

(4.12) 
(4.13) 



In order to solve the eigenequation ( l.lOj ), we observe the highest weight state corresponding to M down-spins and 
N — M up-spins, 



M 



(4.14) 



where 



'++■■■+ 



stands for the state of N up-spins, s m flips a up-spin to a down-spin at the m^th position. The state 



(4.14 ) is checked to be a common eigenstate of both operators Tj and Tj if ip takes the following form 

M 

f mi m 2 - mM = T ^ II F ( (°^)/*> m /J 

(7 fj,— 1 



(4.15) 



G 



where 



r(o-jQ-A) _ [i(d-\) i+1 - (crA)j - c] T(a N aX) 
r(erA) ~ [t(erA)<+i - (erA)< + c] ' r(crA) 



-1. 



and 



F(X; m) = 



i(kj - A) - c/2 
i(fc m - A) + c/2 11 - A) + c/2' 



n 



(4.16) 



(4.17) 



where a stands for the elements of the Weyl group of the Lie algebra Bm] = 1, 2, • • • , M) and <jm stands for its 
basic elements. The parameter A M (/x = 1, 2, • • • , M) for the spin wave function is supposed to satisfy 



M 

n 



Au - A„ + ic A„ + X v + ic 



A„ - Ay - ic A„ + Ay - ic 



n 



fcj — ic/2 A M + fcj — ic/2 



(4.18) 



1 A M — fcj + ic/2 A M + fcj + ic/2 
This is a generalization of Bethe-Yang ansatz The eigenvalues corresponding to the operators Tj and T^ 1 are 



M 

n 



A^ — fcj — ic/2 

1 1 A w - kj + ic/ 

w= i j 



a/ 



_x tt Ay + fcj — ic/2 
2' Zj = 11 Ay + fc,- +ic/2 



li - Ay + kj + ic/1 



respectively. Hence we get the coupled Bethe ansatz equations for the variables A M , (/i = 1,2, 
quasi- momenta Zjfj" 1 — /?, explicitly 



K j ikj -2ikjL 
. . c 



ic/2 Ay + kj — ic/2 



M \ 

11 Ay - kj + ic/2 Ay + kj + ic/2 



Taking the logarithm of (4.20) and (4.18) we obtain the secular equations 



M r 



*-z'<-z-- 1 <&-z£ 



tan -i(^zA ) + tan -i(4l±^ 



c/2 



c/2 



(4.19) 
,M).: and the 

(4.20) 



(4.21) 



E 



V c/2 



V c/2 



M 



tan- 1 ( ) + tan - x ( K+± • 
c c 



(4.22) 



where both and J p are integers, that play the roles of quantum numbers. These equations determine the spectrum 
{kj} of the system with total spin S z = — M and the total energy E — J2f=i tf- 



V. CONCLUSIONS AND REMARKS 



In the previous sections, we solved the problem of N fermions with <5-function interaction in a one dimensional 
potential well of finite depth. In our discussion, we considered the special case in which a single particle tunnels out 
of the potential well. In this case, we showed that there exists a Bethe ansatz like exact solution for the problem. 
We also obtained the secular equations which determine the spectrum of TV-body system. Differing from the usual 
(5-function problem, there is a contribution to the secular equations by the reflection of the particles at the ends of the 
potential well. This contribution affects both th e cha rge and spin excitations. A new feature of present model is the 
additional term proportional to sin~ 1 (kj /Vq) in ( 4.21 ). As a phase shift of one particle at the barrier of potential well 
of finite depth, this term contributes to charge excitation only. The phase shift vanishes when Vq — > oo, which agrees 
with the result for the potential well of infinite depth. In the thermodynamic limit [L — > oo, iV — > ooandM — > oo 
but N/L — M/L keeps a finite value), equations (4.21) and (4.22) will become integral equations. Certainly, one can 
discuss the ground state properties and its thermodynamics in the spirit of |1^| . 

We known that the N particles were thought of as situated on a circle in the case of periodic boundary conditions. 
For the problem of N electrons moving in a quantum wire of finite length, the periodic boundary condition is no 
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longer available. Then we may consider the boundary as a potential well. If all the TV electrons are confined in an 
interval, the potential well is supposed to have infinite depth. Strictly speaking, electrons are able to leave the ends 
of the wire due to the tunneling effect. Then the more appropriate boundary condition is a potential well of finite 
depth. Different from the case of infinite depth in which one needs only to consider the wave function in the potential 
well, we should take into account the possibility of a non-vanishing wave function outside of the well in the case of 
finite depth. The single particle tunneling process is the simplest non-trivial one for the problem of the potential well 
of finite depth. The result of present paper is the first step toward the exact solution of this problem. 
This work is supported by NNSF of China and NSF of Zhejiang province. 
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